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In this paper, we study the bulk motion of a classical extended charge in flat spacetime. A 
formalism developed by W. G. Dixon is used to determine how the details of such a particle's 
internal structure influence its equations of motion. We place essentially no restrictions (other 
than boundedness) on the shape of the charge, and allow for inhomogeneity, internal currents, 
elasticity, and spin. Even if the angular momentum remains small, many such systems are found 
to be affected by large self-interaction effects beyond the standard Lorentz-Dirac force. These are 
particularly significant if the particle's charge density fails to be much greater than its 3-current 
density (or vice versa) in the center-of-mass frame. Additional terms also arise in the equations of 
motion if the dipole moment is too large, and when the 'center-of-electromagnetic mass' is far from 
the 'center-of-bare mass' (roughly speaking). These conditions are often quite restrictive. General 
equations of motion were also derived under the assumption that the particle can only interact with 
the radiative component of its self-field. These are much simpler than the equations derived using 
the full retarded self-field; as are the conditions required to recover the Lorentz-Dirac equation. 

PACS numbers: 03.50.De, 04.25.-g 



Originally motivated by the discovery of the electron, the behavior of small (classical) electric charges has been 
studied in various contexts for over a century. The first results were obtained by Abraham Q for a non-relativistic 
non-spinning rigid sphere. This calculation was later repeated by Schott and Lorentz within special relativity, 
and by Crowley and Nodvik [4j in general (background) spacetimes. Similar results have since been obtained by a 
number of authors for more general charge distributions Ifj, kg, M , l8l l9l llfl. lll| in flat spacetime. These results have 
very recently been extended to curved backgrounds as well |39j. Still, these derivations did not allow for significant 
elasticity, charge-current coupling, and/or rotation. 

Detailed reviews of different aspects of the self-force problem (in electromagnetism as well as scalar field theory and 
general relativity) have been given by Poisson 01 > Havas [3, and Spohn To summarize, though, a common 

theme throughout all of these works has been that the equations of motion describing sufficiently small particles were 
found to be independent of the details of their internal structure to a considerable degree of precision. Only a few 
parameters such as the rest mass and total charge entered into these equations. Ignoring the effects of spin and 
internal currents, the apparently universal correction to the Lorentz force law is given by the well-known Lorentz- 
Dirac self- force. If a particle with charge q has a center-of-mass position z a (s) (with s being a proper time), then this 
self- force is equal to (using a metric with signature —2, and units in which c = 1) 



Although this is only an approximation to the full self-force for any realistic (extended) charge, it is natural to 
introduce a new concept into the theory for which it is exact - that of a point particle. The immediate problem with 
this is of course that the self-field diverges at the location of any point-like source, which would appear to imply that 
its self-force and self-energy are not well-defined. Dirac 0] removed this problem by noting that self-forces acting 
on a structureless particle should only arise as a reaction to emitted radiation. Letting F£^ { , + s denote the advanced 

self-field and F^ l{ ,\ the retarded one, we can define 'radiative' and 'singular' portions of F^ 1{ ,_-. = F^ l{(S s + F^ U(R) : 



I. INTRODUCTION 
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As the names imply, the singular field contains the entire divergent part of the retarded field. It is also derived from 
a time-symmetric Green function, so it would not be expected to contain any radiation. The self-force on a point 
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particle should therefore be determined entirely by the radiative self-field. Combining it with the Lorentz force law 
immediately recovers (ll.lf) 12] . 

This prescription is much simpler than any direct derivation of the equations of motion for a finite charge, and for 
this reason, it has been generalized to work in curved spacetime, as well as for scalar and gravitational self-forces 
[r7ll34| (although this is not the only way of 'renormalizing' point particle self- fields [lill^j). As before, the relevance 
of these extensions to realistic extended bodies has also been established in certain special cases. In linearized gravity, 
for example, a small nearly-Schwarzchild black hole has been found to obey the same equations of motion as a point 
particle |17ll25l |. More generally, it has been shown that a nonspinning body's internal structure is irrelevant to very 
high order within Post-Newtonian theory 0| . Still, there remain questions of exactly how universal these results are. 
Will a spherical neutron star in near-equilibrium fall into a supermassive black hole in the same way as another one 
that is spinning rapidly and experiencing internal oscillations |38f : or one having a mountainous (solid) surface |37|? 
Such systems could be important sources for the upcoming Laser Interferometer Space Antenna (LISA). 

Rather than addressing such questions directly, we have chosen to study the motions of charged bodies in flat 
spacetime (in part) as a model problem. The methods used here were specifically chosen so that very few conceptual 
changes would be required to consider particles (charged or not) moving in a fully dynamic spacetime. This has 
led to some additional complexity not strictly necessary to solve the problem at hand, although the majority of the 
complicating issues have been placed in the the appendix. 

Our derivation is based on W. G. Dixon's multipole formalism pol l2ll I22I I23L |24| . This gives a relatively simple, 
unified, and rigorous way of understanding the motions of arbitrarily structured bodies in both electrodynamics and 
general relativity. Despite the fact that Dixon's theory decomposes the source functions into multipole moments, 
we never ignore any of them. By using generating functions, the entire infinite set is retained throughout all of our 
calculations. This only appears to be possible in this formalism. The laws of motion that are used are therefore exact 
(despite being ordinary differential equations) . All of the approximations that we make are only used to compute the 
self-field. 

Taking a hybrid point of view where an extended body can only interact with its radiative self-field, it is found that 
the Lorentz-Dirac equation follows for a wide range of nonspinning charges with small dipolc moments. If the body 
instead interacts with its full retarded self-field, this result no longer holds. In this case, the equations of motion are 
drastically different if the charge and (3-) current densities have remotely similar magnitudes in the center-of-mass 
frame. Although it is impractical to define exactly what this means at this point, it will be shown that one of these 
quantities must be at least 'second order' compared to the other in order for these extra terms to vanish. This is 
because these cases usually allow the self- fields to do significant amounts of internal work (e.g. Ohmic heating). Even 
when the charge-current coupling in negligible, there are still extra complications in the retarded case. The conditions 
required to generically exclude these and other complicating effects are derived, and turn out to be surprisingly 
restrictive. 

Sec. [n] reviews the various steps involved in calculating the motion of matter interacting with an electromagnetic 
field. Sec. IIIII then summarizes the appropriate definitions of the center-of-mass and its laws of motion as obtained 
from Dixon's formalism. It also decomposes the 4-current in a particular way that happens to be convenient in this 
framework. Although this reduction is not strictly required for the current problem, it is adopted throughout on the 
grounds that it would be essential in curved spacetime. It is derived in detail in the appendix, which also contains an 
in-depth review of of Dixon's ideas. 

With these basic ideas in place, Sec. IIVI goes on to derive expansions for the advanced and retarded self-fields of 
an arbitrarily-structured charge using a slow motion approximation. Sec. Ivlthcn combines these results with those 
from Sec. IIIII to find general expressions for the self- force and self-torque. Finally, Sec. IVII examines the equations of 
motion for certain simpler classes of charges, and derives some conditions under which the Lorentz-Dirac equation is 
applicable. 

We use units in which c = 1 throughout. In order to facilitate a more direct comparison to Dixon's papers, the 
metric is chosen to have signature —2 (although the rest of the notation used here frequently differs from Dixon's). 
We also assume that the spacetime is flat, and adopt Minkowski coordinates for simplicity. Latin indices refer to these 
coordinates, while Greek ones are triad labels running from 1 to 3. 

II. THE PROBLEM OF MOTION 

In studying the dynamics of any system in a classical field theory, one has to specify 'laws of motion' for both the 
field and matter variables. In our case, the only field is of course the electromagnetic one, F ah — . As usual, this 
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is governed by Maxwell's equations 

d b F ab = J a , (2.1) 

glapbc] = q _ ( 2 2 ) 

We assume that the matter in our problem is completely described by its stress-energy tensor T ab and 4-current 
vector J°. Taking the divergence of (|2.1|1 immediately gives our first constraint on these quantities: 

d a J a = . (2.3) 

This equation acts as one of the laws of motion for the matter fields. The other is derived from the requirement that 
d a (T ab + T^A — 0, where is the stress-energy tensor of the electromagnetic field. Combining the standard form 
of |3] with Maxwell's equations then shows that the matter moves according to 

d b T ab = ~F ab J b . (2.4) 

(|2.1() - (|2.4|l are essentially the entire content of classical continuum mechanics in flat spacetime. Of course, different 
types of matter do not all move in the same way, so these equations by themselves are not sufficient to determine T ab 
and J a for all time (even if the F ab were given). One also needs to specify something analogous to a (generalized) 
equation of state, which can take on a rather unwieldy form. 

This sort of procedure is the standard one in continuum mechanics. Unfortunately, the resulting nonlinear partial 
differential equations are notoriously difficult to solve. Such a detailed description of the system should not really be 
required, however, for problems where we are only interested in the body's bulk motion. In these cases, a representative 
world line could be defined inside the (convex hull of the) spacelike-compact support of T ab . Given that this worldline 
can always be parametrized by a single quantity, its tangent vector might be expected to satisfy an ordinary differential 
equation - at least when using certain approximations. Solving such an equation would clearly be much more 
straightforward than the partial differential equations that we started with. 

This sort of simplification is one of the main motivations behind the many (source) multipole formalisms in the 
literature 0, Isflj . In these, one first fixes some particular reference frame which has, among other properties, a 
preferred time parameter s. The quantity being expanded - say J a (x) - can then be written in terms of an infinite 
set of tensors depending only on s: Q a (s), Q ab (s), . . . The reverse is also true. Given J a , there is a well-defined way 
to compute any moment. The set {Q "} is therefore completely equivalent to J a . This implies that the conservation 
equation 1)2. 3|) may be used to find restrictions on the individual moments. 

Such restrictions depend on the precise definitions that are being adopted, but can usually be divided into two 
general classes. The first of these consists of purely algebraic equations imposed at a fixed value of s. We call these 
the constraints. There are also a number of evolution equations which usually take the form of ordinary differential 
equations. Multipole expansions can therefore be used to convert (|2.3|l and l|2.4|l into a number of algebraic and 
ordinary differential equations (without any approximation). 

This does not actually simplify things as much as it at first might appear. The reason is that almost all definitions 
for the source multipolcs will lead to an infinite number of (coupled!) evolution equations. This is often dealt with 
in practice by assuming that all moments above a particular order are irrelevant, which leaves one with only a finite 
number of evolution equations. There are, however, interesting questions that require knowing the higher moments. 
J a (or T ab ) for example, cannot be reconstructed without them. This means that the self-field cannot be calculated in 
the near zone from only the first few moments. Although the self-force and self-torque can be found in certain cases by 
examining energy and momentum fluxes in the far zone 

[3 13 13 El (which can be adequately approximated using 

only a finite number of moments), this is considerably less accurate than integrating the force density throughout the 
charge's interior. 

For these reasons and others, it is desirable to define a set of multipoles that do not require any cutoff. Remarkably, 
such a set exists [3] , and is essentially unique [3] • Without any approximation, moments for both J a and T ab can be 
defined which satisfy a finite number of evolution equations. There remain (uncoupled) constraint equations for each 
moment, although these are easily solved. We adopt this formalism due to Dixon for the remainder of this paper. Its 
net effect is to allow us to relate the motion to the fields in a more rigorous way than is usually done (short of directly 
solving 1)2. 3|) and l|2.4|l ). It does not, however, have anything to say about the fields themselves. We therefore obtain 
the self-field in a standard way, and then use Dixon's equations to find how the matter moves in response to it. 



III. LAWS OF MOTION 

As noted, we use Dixon's method [20l l2ll [3 13 I3| to decompose J a and T ab into multipole moments. Each of 
these sets is designed to describe as simply as possible all possible forms of J a and T ab satisfying their respective 
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conservation equations. We first assume that these matter fields are at least piecewise continuous, and have (identical) 
supports W . Any intersection of a spacelikc hypcrsurface with this worldtube is assumed to be compact. 

Now choose a timelike worldline Z C W, and a timelike unit vector field n a (s) defined on Z. It is assumed that 
this is always possible in any physically interesting system. Z is then parameterized by the coordinate function z a (s), 
and the tangent vector to it is denoted by v a (s) :— dz a /ds =: z a (s). v a need not be equal to n a , although it will be 
convenient to normalize s such that n a v a = 1 (so v a v a 7^ 1 in general). v a is called the kinematical velocity, while n a 
is the dynamical velocity. The set {Z, n a } then defines a reference frame for the definition of the multipole moments. 
At this point, it should be thought of as arbitrary, although physical conditions will later be given that pick out a 
unique 'center-of-mass' frame. 

A collection of spacelike hyperplanes {£} can easily be constructed from Z and n a . Each £(s) is to pass through 
z a (s), and be (everywhere) orthogonal to n a (s). Assume that any point in W is contained in exactly one of these 
planes. Unless n a is a constant, it is clear that this property cannot be true throughout the entire spacetime. For 
x G S(s) n W, we must therefore have that ma,x\h a (s)(x — z ( s )) a \ < 1 (among other conditions), which gives a weak 
restriction on the body's maximum size. It is not really important, though, as any reasonable type of matter would 
be ripped apart long before this condition was violated. 

The main results that we need from Dixon's theory at this point are his definitions of the linear and angular 
momenta. These disagree with the usual ones when either J a or F ab are nonzero, although it is still convenient to 
label them by p a (s) and S ab (s) respectively. Unless otherwise noted, the words 'linear and angular momenta' will 



always refer to the quantities |20| |21| 

P a (s) = / dE 



b 

E( S ) 

S ab (s) = 2 / d£ c 



T b 



T ab + J"r c / du F ac (z(s) + ur) 







[a T b]c + J c r d l duur [a F b]d (z{s) 







(3.1) 
(3.2) 



where r a := x a — z a (s), and u is just a dummy parameter used to integrate along the line segment connecting z a (s) 
to x a = z a {s) +r a . 

Detailed motivations for these definitions can be found in as well as the appendix. In short, though, it 

can be shown that the given quantities are uniquely determined by demanding that stress-energy conservation directly 
affect only the first two moments of T ab (once the concept of a moment has been defined in a reasonable way). If p a 
and S ab are known in some time interval and satisfy the appropriate evolution equations, the class of all stress-energy 
tensors with these moments can be constructed without having to solve any differential equations. Each of these will 
exactly satisfy l|2.4l) . 

This property implies that evolution equations for the quadrupole and higher moments of the stress-energy tensor 
are nearly unconstrained. They can be thought of as the 'equation of state' of the material under consideration. This 
type of independence of the higher moments from the conservation laws also occurs in Newtonian theory 0| > an d 
there are considerable advantages in preserving as much of that structure as possible in the relativistic regime. In 
particular, there is no need to discard multipole moments above a certain order. The choices l|3.1[) and l|3.2|l allow us 
to retain many of the conveniences of a multipole formalism without its classic limitations. 

The same definitions for the momenta can also be motivated by considering charged particles in curved spacetime 
|2lL |24|. There, one can study the conserved quantities associated with Killing vectors in appropriate spacetimes 
(where both the metric and electromagnetic field are assumed to share the same symmetries). Fixing Z and n a (s) 
allows each such quantity to be written as a linear combination of vector and antisymmetric rank 2 tensor fields on 
Z. Crucially, the definitions of these quantities do not depend on the Killing vector under consideration, so we can 
suppose that they are meaningful even in the absence of any symmetries. If the metric is taken to be flat, these 
objects reduce to l|3.1|) and l|3. 2f> 21] . They are therefore the natural limits of what would generally be referred to as 
the linear and angular momenta of symmetric spacetimes. 

It is now convenient to define a coordinate system on W that is more closely adapted to the system than the 
Minkowski coordinates used so far. First choose an orthonormal tetrad {n a (s), e° (s)} (a = 1, 2, 3) along Z. Requiring 
that each of these vectors remain orthonormal to the others implies that 

e a a = -n a n b e b a , (3.3) 

which is essentially Fermi- Walker transport. A spatial (rotation) term may also be added to this, although our calcu- 
lations would then become considerably more tedious. On its own, the choice of tetrad has no physical significance, 
so we choose the simplest case. 

Since it was assumed that {£} foliates W, any point x € S(s) fl W may now be uniquely written in terms of s and 
a 'triad radius' r a 

x a = e a a (s)r a + z a (s) . (3.4) 
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Varying r a with a fixed value of s clearly generates S(s). Also, the Jacobian of the coordinate transformation 
(x a ) — > (r Q , s) is equal to the lapse iV of the foliation, 



iV(a;) 



n a (s)r° , 
n a (s)e a a (s)r a 



(3.5) 



We will be extensively transforming back and forth between these two coordinate systems, so it is convenient to 
abuse the notation somewhat by writing f(r a ,s) — f (x a ) for any function /. The intended dependencies should 
always be clear from the context. It is also useful to denote quantities such as n a e^ by h a . Note that in this notation, 
n a = n Q e£ ^ dha/ds. 

Using these conventions, it is natural to split J a into the charge and 3-current densities seen by an observer at 
r = (the 'center-of-mass observer') 



J a = pn a + e a J a . (3.6) 
It is then shown in the appendix that there exist 'potentials' ip and H a which generate J° through the equations 

P = 9 Q (r» , (3.7) 
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N~ 



H a + v f3 d f3 (r a ip)-r a p 



(3.8) 



Denoting the total charge by q and letting |r| 



-r a r a > 0, ip(r a ,s) was found to be continuous, and equal to 



q/A7r\r\ 3 outside W. Similarly, H a (r^,s) is given by 1A40|) outside W, and is piecewise continuous in r' 3 . H a also 
satisfies d a H a = 0. These properties guarantee that J a has support W, and is everywhere continuous. A direct 
calculation also shows that 



d a J a 



dp3 







= 



(3.9) 
(3.10) 



as required by l|2.3|l . 

Any physically reasonable current vector can now be constructed by choosing potentials satisfying these rules. A 
physical interpretation of one's choice is then given by substitution into (|3.7|l and l|3.8|l . For example, a uniform 
spherical charge distribution with time-varying radius D(s) is described by (assuming n a = v a ) 



<p(r, s) 
H a (r, s) 



Q(D(s 



\r\) + (^f) 3 &(\r\-D(s)) 



4ttD 3 (s) 
, 

where 8(-) is the Heaviside step function. H3.7|l and Ij3.8|l then show that the tetrad components of J a are 



P(r, s) 



3 a (r,s) = 



em 



3qr a 



D(s) 

4TrN(r,s)D 3 (s) \D(s) 



Q(D(s)-\r\) 



(3.11) 
(3.12) 

(3.13) 
(3.14) 



It is also clearly possible to calculate ip and H a from any given J a . This requires inverting l|3.7|l . which acts 
as a partial differential equation for ip on each time slice. The solution to this equation would usually have to be 
obtained numerically, which is clearly inconvenient. Largely for this reason, we shall consider tp and H a to be the 
given quantities for the remainder of this paper. J a can be derived from them using operations no more complicated 
than differentiation. 

Another reason for this unconventional choice is that ip and H a contain all of the multipole moments of J a in a 
natural way. As shown in the appendix, they are closely related to the Fourier transform of a generating function 
for these moments. An arbitrary current moment can be obtained essentially by differentiating the inverse Fourier 
transforms of the potentials a suitable number of times. For example, (|A1|) . Ij A24|l . I|A30(1 . (|A3ip . and (|A38|I can be 
used to show that the dipole moment has the general form 



= d 



2n^ a e%PN ( ip(r, s) 



4vr I r- 1 3 



■e a a e b H^(r,s) 



(3.15) 
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where H af3 = is defined by 



Given (IA40I). we let 



H a = d B H ap 



H 



2ttH 3 



outside W. 

For the example given in (|3.11|) and H3.12[) . the dipole moment is equal to 



Q ab = -qD z {s)n^h 



1 



K6] 



(3.16) 



(3.17) 



(3.18) 



This might have been expected to vanish in spherical symmetry, although it should be noted that it is only nonzero 
when sphericity is being defined in an accelerated reference frame. 

These constructions can now be used to simplify the definitions of p a and S ab . Letting f a := ur a , the electromagnetic 
term in 1|3.1|) is equivalent to 



J d 3 fe^f F ab (f,s) dittt -4 p(*7u. s) . 



But H3.7fl shows that 



duu 4 p(r/u,s) 



ou 



tp(r, s) 



u 3 I <p(r/u, s) 



47r|r/up 



47r r ' 



so (|3.1(l and (|3.2(l can be rewritten as 

p a = I d 3 r 

S ab = 2 / d-V 



n b - \<p — 



r~,elF ac 



47r|r| 3 J 7 c 



(3.19) 

(3.20) 
(3.21) 

(3.22) 
(3.23) 



We now need evolution equations for these quantities, which are most easily obtained by direct differentiation. For 
any function I b (x), relating the ^-coordinates of (r a , s + ds) to those of (r a , s) shows that 



d 

dl 



dE b I b {x) = / d 3 r h b I b + {v c - n c h a r a ) d c (n b I b ) 



If I b vanishes outside of some finite radius, this expression simplifies to 



_d 

ds 



/ d£ b / b (x) = / d 3 rNd b I b 



(|2.4|) . I|3.22ll . and (|3.23|) can now be used to show that 



and 



d 3 r ^ NF ab J b 



4^|r|3 J 7 c 



r^e7F° 



(3.24) 



(3.25) 



(3.26) 



S ab = -2 1 d 3 r I Nr a e\?F b ^J c + v^ a T b ^ c n c + 



r^eJr a e^F b ^ c 



47r | r |3 



= 2p^ a v b] - 2 y d 3 r J iVr^'^^J, + -|- 



47Tk|3 ' C 



r^elv^F^ 



(3.27) 
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If the momenta had been defined in the usual way, the (ip — g/47r|r| 3 ) terms would be absent from these expressions. 
The extra complication in the present case derives from the electromagnetic couplings in (|3.1[l and (|3.2[) . 

Unsurprisingly, H3.26|l and (|3.27() can be directly related to the monopole and dipole moments of the force density 
F ah Jb. Denoting such moments by \E' a (s) and ^ ab (s) respectively, it is possible to prove l|A56|) and (|A57|) . It is then 
natural to refer to — vE» a as the net force, and —2^^ as the net torque acting on the body. Viewing these quantities 
as force moments allows one to derive l|A61l) and l|A62|l . But these results are no different than (|3.26|1 and (|3.27|) . If 
desired, the reader may therefore view *° and *[ ab l to be defined by (|A"56)) and (|A"57|) . 

In cases where F ab varies slowly in the center-of-mass frame (both spatially and temporally), the expressions for 
the force and torque can be expanded in Taylor series involving successively higher multipole moments of the current 
density (denoted by Q "). Deriving such equations would be awkward using the methods introduced in this section, 
so we simply refer to <|A63|1 and (|A64|) . Despite the peculiar definitions of the current moments being used, these 
expansions are exactly what one would expect out of a multipole formalism, and can therefore be considered a check 
that Dixon's definitions are reasonable. 

In most cases where (|A63|) and l|A64|l are any simpler than the exact expressions for the force and torque, only 
the monopole and dipole moments will be significant. These moments can be computed from (|AI3|) and (|3.15|1 
respectively. In the limit that the particle is vastly smaller than any of the field's length scales, only the monopole 
moment will enter the equations of motion. In this case, the force reduces to the standard Lorentz expression and the 
torque vanishes, as expected. 

The final ingredients required to complete this formalism are unique prescriptions for Z and n a . Simply knowing 
the linear and angular momenta at any point in time does not necessarily determine the body's location in any useful 
way. The problem is compounded by the fact that these quantities are strongly dependent on the choice of reference 
frame itself. Indeed, without knowing where W is, it is essentially impossible to specify J a in any meaningful way. 

These problems can be removed by choosing Z and n a appropriately, and then assuming that the resulting worldline 
provides a reasonable representation of the body's 'average' position. Following [21H23I I24I |2S | . we first assume that 
for any point z S W, there exists a unique future-directed timelike unit vector n a (z) such that 

p a {z; n) = M(z; n)n a (z) , (3.28) 

for some positive scalar M . Here, we have temporarily changed the dependencies of p a and n a for clarity. It is seen 
from l|3.I[) that p a depends nontrivially on both the base point z, and on n", which defines the surface of integration. 
(|3.28l) is therefore a highly implicit definition of n a . 

In any case, another condition must be also be given to fix z. We want this to lie on a 'center-of-mass line' in some 
sense, so it would be reasonable to expect the 'mass dipole moment' defined with respect to it to vanish: 

n b n c t abc = , (3.29) 

where t abc is the full dipole moment of the stress-energy tensor. This is defined by 1A47(I . so H3.29[l is equivalent to 

n a (z)S ab (z;n) = . (3.30) 

The integral form of l|3.30|l reduces to the standard center-of-mass condition when F ab = and n a = v a . It also 
allows us to write the angular momentum in terms of a single 3-vector S a 

gab = ^bed^g^ (3 31 ) 

Of course, both of these properties would also have been satisfied by instead requiring v a S ab — 0. We reject this 
choice due to the fact that it leads to nonzero accelerations even when F ab = [13 • Replacing v a by n a avoids 
this peculiar behavior, which we consider to be an important requirement for anything deserving to be called a 
center-of-mass line. It also seems more natural to define the mass dipole moment in terms of n a rather than v a . 

For the remainder of this paper, we assume that l|3.28|l and l|3.30|l are always satisfied, and call the resulting Z and 
n a (s) the center of mass frame [2l|, [24], [28| . It is not obvious that solutions to these highly implicit equations exist, 
although existence and uniqueness have been proven in the gravitational case |29| . There, it is also true that (given 
certain reasonable conditions) Z is necessarily a timelike worldline inside the convex hull of W. We assume that the 
same results extend to electrodynamics. 

The uniqueness of these definitions is actually not very critical for our purposes. The important point is that 
a solution with the given properties can presumably be chosen for a sufficiently large class of systems. While it is 
clearly very difficult to find the center-of-mass directly from T ab , J a , and F ab , it is relatively straightforward to simply 
construct sets of moments which automatically incorporate i|3.28[l and (|3.3U|) . 1)3. 7|) and 13.8fl show, for example, that 
these definitions do not have any effect on our ability to construct arbitrary current vectors adapted to them. There 
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is nothing preventing the moment potentials ip and H a from being appropriately centered around r a = 0. Although 
it is not obvious that this can also be done for the stress-energy tensor, we conjecture that it can. 

Now that H3.28|) and l|3.30JI have been assumed to hold, we need evolution equations for M, n a , and z a . These are 
easily found from (|A56|I and l|A57|) : 

M = -n a ^ a , (3.32) 

Mh a = -h^ b , (3.33) 

M(v a - n a ) = S ab h b - 2^n b , (3.34) 

where h% is the projection operator 5% — n a n b . Note that the last of these equations shows that n a = v a if the spin 
and torque both vanish. 

In general, l|A"57)) and IpO^ - lTOH may be used together to find the motion of the body's center-of-mass in terms 
of \l/ a and \J/[ ah l. Once the field is known, these quantities follow from (|A61|) and (|A62|) . Some recipe for evolving the 
dipole and higher current moments in time - most conveniently expressed in terms of ip and H a - is also required. 
Combining all of these elements together leaves us with a well-defined initial value problem that will determine z a , 
n a , M, S ab , and J a . If the stress-energy tensor is also desired, possible forms of it could in principle be constructed 
from 1|A51|) and (|A54() in the same way that J a was derived from (p and H a . Combining all of these steps would 
completely characterize the system, although we shall omit the last one. The result is still sufficient to answer most 
questions that one would be interested in asking of a nearly isolated particle. 

IV. THE FIELD 

It is clear from the previous section that the equations of motion will easily follow once \I/ a and <]>[ ab ] are known. 
These depend on the field, which we now calculate. For a reasonably isolated body, it is first convenient to split F ab 
into two parts 

F ab = F c ab t + F s f u . (4.1) 

The external field is assumed to be generated by outside sources, while the self-field is entirely due to the charge itself, 
bmce l|MT|) and (|X62)l are linear in F ab , the force and torque can also be split up into 'self and 'external' portions. 
In Lorenz gauge, the vector potential sourced by the particle is given by [8j 

AU(x) = - J dV J a (x')S(a(x, x')) , (4.2) 

where 

a(x,x') = ^{x-x') a {x~x') a , (4.3) 

is Synge's world function [rH l3"H - In realistic systems, this potential will couple to the external one via the outside 
matter fields. These may act as reflectors or dielectrics, or there may be an n-body interaction where the self-fields 
influence the motion of some external charged particles (obviously affecting the fields in W). Although it would be 
reasonable to group together all portions of F ab causally related to our particle in some way as the 'self- field,' this 
would be impossible to do with any generality. Instead, we simply define the F s "j f to be the field derived from 14.21) 
in the usual way (this differs from the point of view taken in e.g. [3f|)- The interactions of the self-field with the 
outside world will all be categorized as 'external' effects that we presume can be accounted for by separate methods. 
If all of the external matter is sufficiently far away from W (and slowly varying), F£ b t will be approximately constant 

within each S(s) n W slice (and from one slice to the 'next'). 5'° xt and ^g"? can therefore be approximated by (|A63|) 
and (|A64|) in these cases. Finding the self-force and self-torque is more complicated. For this, F® b l{ has to be combined 
with the exact integral expressions for the force and torque - l|A61|l and (|A62(1 . The detailed structure of the self-field 
must therefore be known throughout W. 



This easily follows from A* 



self- 



F s ab u (x) = 2d^A b J cli , 

= -2 fd 4 x'S'(a){x~x')^J b ^ . (4.4) 
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Writing x in terms of (r Q , s), and x' in terms of (r' a , t), and defining &(x, x') := da/dr, F™fe becomes 

~N , 





r i d 


j dV < 









-(x-a;') [a J bl 



(4.5) 



t + (> s) represents the advanced time, and r_ the retarded one. These are found by solving a(x,x') — with x and 
r' held fixed. 

Although we only consider the retarded field to be real, the advanced solution is also retained for now. This allows 
us to later construct the radiative self-field, which is considerably simpler than the full retarded field. It would be 
quite convenient if the self-forces generated by the two fields were identical (as Dirac assumed for point particles |12fl'). 
although we will show that this is not true in general. 

Returning to the explicit form for F® b if , splitting up J a according to (|3.6[) shows that 



F£ { (x) = 2 /dV-^-Jp 



(x - x') [a h b ^ - (? + iV -1 n a rJ,J (x - x') [ V] - v [a n b ^ 
(x - a;') [a e bl + N^h^r'^j + n a (x - x')^ a n^ + (Nn^ a + (h ■ v) [a ) e\\ 

+r(x-x') [a e si ■ 



+ p(x - x') [a n b ^ 



(4.6) 



Here, j a (r\r) := dj a (r' ,t)/8t, which differs from our usual convention (e.g. n Q (r) := e"(r)dn a (r)/dr). 
Moving on, l|4.3|l implies that 



&(x, x') = - \N(r', t>°(t) + (ft • u) a (r)) (x - x') a 



and 



a(x,x') = N 2 +v a v a 



(4.7) 



(4.8) 



If we now specify how ip and H a vary in time (which determines p and j a ), we would have all of the ingredients 
necessary to find the body's motion without any approximation. Unfortunately, inserting (|4.6I) into l|A57() . I|A61(I . 
(|A62|I . and (|3.32|l - (|3.34|l leads to a set of delay integro-diffcrential equations for the object's motion. Such a system 
would be extremely difficult to solve (or even interpret) in general, although it could be a useful starting point for 
numerical simulations. It might also be interesting analytically if one were looking for the forces required to make a 
body move in some pre-determined way (e.g. a circular orbit). 

Such questions will not be discussed here. We instead consider the simplest possible approximations that allow us 
to gain insight into a generic class of systems. In particular, it is assumed that all of the quantities in l|4.6|) which 
depend on t± may be written in terms of quantities at s (via Taylor expansion). This requires that nothing very 
drastic happen on timescales less than the body's light-crossing time. This is not as trivial a condition as it might 
seem to be (see e.g. |33|), although we will not attempt to justify it. 

Expressing everything in terms of s rather than t± first requires calculating A± := s — t±. t± was defined by a = 0, 
so A± can be found by Taylor expanding this equation in A±. Assuming that both n a and v a are at least C 3 in time 
for s <E [t±,s], 



<(r±) = e-( S ) + A ± n a (s)h a (s) - ^A 2 ± (n a ( S )h a (s) + n a (s)n a (s)) - (^) 
z a (r±) = z a ( S )-A ±v a ( s ) + iA^( s )-iA 3 ± ^( s ) + ^Ai^(d 2) ) , 



(4.9) 
(4.10) 



where £ ±^ and £± ^ are some numbers between t± and s. 



Then (|3.4|) shows that 

(x - x') a ~ e%(r- r') a + A± ^Nn a + (h ■ v) a ^j + ^A 2 ± (n a h a r' a - Nn a ^j + - 



Alii a . 



(4.11) 
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Everything here has been written in terms of s, and all terms involving quantities such as |n| 3 , |ri| 2 , |n|, \h ■ v\ 2 , 
and \v — h\ have been removed. Since |A±| ~ \r — r'|, these terms can be reasonably neglected if |n|72. 3 <C \h\TZ 2 <C 
\h\lZ <C 1, and \v — h\lZ -C \h-v\ <C \h\lZ, where 1Z := max (|r|, \r'\). In a sense, we are expanding up to second order 
in powers of \h\lZ, and up to first order in \v\. 

The requirement \h\TZ <C 1 must hold for all (r, r') pairs, so it is useful to define the largest possible value of 1Z. 
We call this the body's 'radius' 



D(s) 



max r 

£(s)rw 



(4.12) 



Using it, \h\lZ -C 1 implies \n\D <C 1. This technically restricts the allowable size of the charge, although very few 
reasonable systems would actually be excluded. 

Interpreting the relation satisfied by (h ■ v) a (= v a — n a ) isn't quite as simple. Given 1)3. 34|) . it is roughly equivalent 
to assuming that spin effects are present only to lowest nontrivial order. This is not completely accurate, though, and 
more precise statements will be given in the following section. 

In any case, the inequalities following (|4.11(l will be assumed to hold from now on. Using them to expand a = 0, 
we find that 



R l := |r 



A H 



2v a (r - r') a + A ± N(r, s)N(r', s) + ^A 2 ± n a (r + 2r% + -L A| \n\ 2 



(4.13) 



All but the second term here is already 'small,' but not quite negligible in our approximation. The zeroth order 
expression for A± (= =fR) may therefore be substituted into each of these terms without any overall loss of accuracy. 
The resulting equation is easily solved: 



A± ~ T R 

±R- 1 v a (r -r') a 



1 + \n a (r + r'T - \n a vPv a r' & ± ^Rn a (r + 2r% + ^R 2 n a h 



a + l(ii a (r + r') a y 



(4.14) 



where everything is evaluated at s. 

Although it would still be straightforward at this stage to compute the exact error in l|4.14|l . it is not necessary. 
Dimensional analysis shows that the neglected terms have magnitudes |n|7?. 4 , |ri| 2 7?. 5 , \d(h ■ v) / ds\lZ 2 , \d 2 (h-v)/ds 2 \R? , 
\h ■ v\ 2 lZ, and so on (where each of these can be evaluated anywhere in the interval [t±, s]). 

Continuing to expand quantities appearing in H4.6f) . a long but straightforward calculation shows that (|4.7II and 
(14. 8|) can be approximated by 



a ~ ±R 



1 - X -h a (r + r') a - l(n a (r - r 1 )^ T ^Rn a {r + 2r% + ^i? 2 |rV 



and 



a ~ 1 - h a (r + r') a t Rn a (r + 2r') a + n a h p r a r' + -R 2 \h\ 2 



(4.15) 



(4.16) 



Looking at l|4.6f) . the charge and current densities are the only quantities that have not yet been expanded away 
from t±. It is instructive to temporarily leave them like this, but write out all other terms in field at s. The resulting 
expression is quite long, so we break it up into several smaller pieces by writing F^,^. in the form 

*2Sf(±) = 2 / (p(rV±)/(f) +P(r',T ± )f% -jV,t ± )/$o +3 (r',T±)f$^ . (4.17) 



The definitions of each of these coefficients is obvious from comparison with 1)4. fijl . 

Before computing them, we first simplify the notation by defining a quantity T(s) > such that T^ 1 remains 
(marginally) less than \h\, In] 1 / 2 , \h ■ v^^/D, etc. This is useful because a number of different objects were assumed 
to be negligible in this section. Some of these may be much larger than the others, so using them to writing down 
error estimates would become rather awkward. Introducing T removes this difficulty. 
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Along the same lines, we also define e := D/T <C 1. A simple but exceedingly tedious calculation then shows that 
the coefficients in (I4.17|) are equal to 

f$ ~ -{r-r'Te^ 1 + \h\v - r> - ~ (>(r - r')^) + ^fCrp{r - r% + ±R 2 \h\ 2 

+ ii? 2 (r - r') a e [ Zn b] - ^R 2 n [a n b] (l + ^h (3r - r')p) ± ^R 3 n [a n b] + o (\e{K/D)} "ll ) , (4.18) 



/ab 
0(3) 



/a 



b 



f$ c ±RUr-rTe [ ^ 



1 + h 



±^R 2 (r-r') a e [ ^n b] 



+ (r - r') a e [ ^R (\ + ^(Sr + r')/3^ T ^i? WJ (l + 2n%) - ~rJ a n b] R 3 - (h ■ v)^ a n b ^R 



e(K/D) 



n- 



(r - r'Y^l 1 + \^{r - r') 7 - \{h\r - r') 7 )' + ^h x r y (r - r') x + l -R 2 \n\ 2 
+ ~R 2 n (r - r') Q e[> bl + ]^R 2 np{r ~ r') a e [ > b] + nK b jR 2 Qri 7 (r - r% T ~B\nf 



iT V(r - r% ) + ^R 2 n [a e^ ( 1 + ^n 7 (3r - r% J ± ^R 3 ^ a e 



1 



±R {(r-rTe^ 



e(H/D) 
1 



K 



T ^R 2 (r r'Ye [ S 



+ h^rp + (h%y ± Ji?n /3 (2r + r')/3 - ^# 2 N 
) - (r - rO^n^i* ^1 + in 7 (3r + r') 7 J - n [a e* ] 

- - (n 7 (r + r') 7 ) 2 + n 7 n\ 7 (2r A + r' x ) ± ^i?n 7 (5r + r') 7 - 7^7 # 2 N 2 ± iT V(r - r') 
8 \ / 6 24 



1 + -n 7 (3r-r') 7 



e(K/D) 



(4.19) 



(4.20) 



(4.21) 



A. Approximations 

At this point, it is useful to consider the physical meaning behind our approximations more carefully. Although error 
estimates have been given for for fm-fpuy the important errors are those in F^ lf . Unfortunately, the unbounded 
integrand in (|4.17() makes it difficult to compute these rigorously. We shall simply assume that J a is sufficiently 
homogeneous that dimensional analysis can be used to say, for example, that the absolute value of the error in 
JdV p{r',T±)ffi } /R 3 is less than about 



e 3 max \p\D 



(4.22) 



With this established, we can also use the assumptions outlined in the previous section to provide bounds on the 
force and torque. Combining |n| < T _1 with (|3.33|) . 

< MT- 1 . (4.23) 

Since T 3> D, this means that the particle cannot be accelerated up to an appreciable fraction of the speed of light 
within a light-crossing time - a very reasonable restriction. 
Differentiating (|3.33|l . 

(4.24) 

(4.25) 
(4.26) 



Mii a = - 2Mn a 

Generically, this implies that |fi| < T~ 2 can be ensured by assuming 



M 



< MT~ 



* Q < MT' 
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Differentiating (|3.33|1 a second time and using (|3.32|) . we also have that 

M < AIT- 2 , 
* Q < MT~ 3 , 
n a *° < MT~ 2 . 



(4.27) 
(4.28) 
(4.29) 



if 



Similar restrictions on the spin and torque can be found by using (|3.34() . In particular, |u| will remain less than e 2 



S' 



a/3 I 



\S a \ < eMD 



< e z M 



Differentiating (|3.34|l . \n — v\D <C e 2 is implied by 



< eM , 



* [afc] e> b < eMD 



(4.30) 
(4.31) 

(4.32) 
(4.33) 



Although weaker restrictions than these can be adopted in special cases, we assume for simplicity that they always 
hold. More concisely, we are restricting ourselves to systems where the magnitude of each tetrad component of the 
(full self + external) force is bounded by about MT~ 1 1 and the magnitude of each torque component is no larger than 
e 2 M. Bounds on the s-derivatives of these quantities are suppressed by appropriate factors of T. These conditions 
basically mean that the body can only change significantly over timescales larger than T. I|4.30|l . for example, implies 
that this timescale sets a lower bound on the charge's rotational period. 

These interpretations should not be surprising. Unfortunately, though, they are quite not as simple as they appear. 
Combining them with l|A61|) and l|A62(l can lead to less obvious bounds on the structure of J° itself. One of these 
affects the size of the particle's dipole moment. If the self-torque is small compared to the external one, then the 
magnitudes of 



I 7" r cxt ) ll aCM c^ cx t, 7-^x1;) e a e 6 n c Ud^ -T f 



bd 
ext 



(4.34) 



should all be less than about e 2 M. If q is negligible, these conditions ensure that both the force and torque are 
sufficiently small. 

In most cases of interest, though, q ^ 0. Bounding the quantities in (|4.34l) then guarantees only that the torque 
is acceptable. The (full) force will often be dominated by the monopole component of ^xti in which case (applying 
(l4"2H) and JOd}) 



MT- 1 > |ge> 6 F e 



ab I 
ext i 



(4.35) 
(4.36) 



Since vp is small, this second bound allows the magnetic field to be quite large. Despite this, we choose to assume 

(purely for simplicity) that both the electric and magnetic fields are bounded by M(\q\T) .If the fields are in fact 
as large as these equations allow, then the previous bounds on the dipole moment imply that 



Q a0 \<e\q\D , 



(4.37) 



along with a similar restriction on \Q a e^n^. This implies that the 'center-of-charge' must be very close to the 
center-of-mass, which considerably reduces the number of current distributions that we can allow. Note, however, 
that these assumptions could be relaxed somewhat if extra conditions are placed on the sizes of the field components. 

To find some other implications of our approximations, we now study a particular example in detail before moving 
on to computing general self- forces. 
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V. SELF-FORCES 
A. Stationary Case 

It was shown in the previous section that (|4. 17|) - (|4.21|> approximate the self-field due to a charge with an (almost) 
arbitrarily accelerating center-of-mass line. These equations become exact when n a — and n a = v a . Temporarily 
assuming that these conditions are true, 



(r-rr4Vl(p(r',r ± )T#',T ± )) - (r - r') a ^4 (/(r', r±) 



±Rf{r',T ± )) +n^e^R 2 'f{r' 1 T ± ) 



(5.1) 



where t± = s ± R. 

One might at first think that the self-forces generated by this field would be identically zero, but this is not generally 
correct. Anything emitting a focused beam of radiation, for example, will experience some recoil. It will also lose a 
bit of mass over time, and can even start rotating if the the beam is offset from the source's center-of-mass. These 
can all be thought of as self-force and self-torque effects, and clearly do not require v a ^ 0. 

Since (|5.1|l is exact, one could compute self-forces and self-torques for a variety of systems without any approxima- 
tion. If these were nonzero, the assumptions leading to (|5.1() would not be maintained unless there were also external 
forces and torques present that exactly balanced them. Such a situation would be rather artificial, so we instead 
apply the approximations of Sec. IIVI and use the resulting self-force and self-torque to gain some intuition into more 
general cases. 

Assuming that \<p\ 3> \(p\D 3> |</3|-D 2 , except perhaps at isolated points, one can show that 

*£if(±) - / d 3 r clV ig I R 2 p(r, s)H a (r', s) - 2(r - r)«^(r, a) Qr'^(r', e) T RH^(r', s)j 



tp(r, s) 



4tt r • 



(r-r') a r p -5 a p {r-r') 
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H^r',s)\, (5.2) 



to first order in (p and H a . This illustrates that \E' a is bounded by ipHD 5 (meaning max \<p\ max \H\D 5 ), ipHD 5 , or 
HHD 5 (whichever is largest). If H a is very small compared to ip, then the lowest order contributions to the self- force 
can be shown to be of order ipipD 6 . 

It is then clear that unless the system is particularly symmetric, self-forces with these magnitudes are unavoidable. 
They would considerably complicate the equations of motion when considering nonzero center-of-mass accelerations, 
so it would be convenient to ignore them. This is self-consistent with the approximations already in place if the 
charge's 'elasticity' is bounded by 

ipD < eV , (5.3) 
CpD 2 < eV , (5.4) 
HD < <?H . (5.5) 

If these relations hold, the self- force becomes negligible whenever n a = and n a = v a ; i.e. a charge at rest will 
remain at rest unless acted on by an external field. For this reason, all of the calculations that follow will assume 
(|5.3|) and l|5.4|l to be true. I|5.5|l is extremely simple to relax, so it will be eventually be replaced with HD < e 2 H. 
The interested reader could easily weaken these conditions even further, although we regard the extra complication 
to be unnecessary for the present purposes. 

In writing down i|5.1f) . both h a = and n a = v a were required. This second condition was adopted here for 
simplicity, although it is by no means 'natural.' Given 13. 3411 . it can only hold if V&^rib = 0. But combining (|5.3[l - 
(|5.5() with (|A62|) and (|5.1(l shows that this portion of the torque is usually nonzero even if (p = H a = 0. Generically, 
one finds that ty<- ab >rib will only be small if either ip -C H or H <C <p>. n a ~ v a is not a particularly important 
requirement, so we choose not to impose these restrictions on the current structure. In general, then, n a ^= v a even 
when n a = (and tp can be of order H). 

Of course, \n — v\ — \v\ still can't exceed e 2 . As was already shown, this remains true if each tetrad component 
of vj/[ af> ] is bounded by e 2 M (and |5| < eMD). We now strengthen this assumption slightly, and require the self and 
external components of the torque to be individually less than e 2 M. There are then terms in the self-torque which 
violate this unless more restrictions are placed on the particle's current structure. 
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Using l|5.1[) and (|A62|I again, one can see that the tetrad components of will be of order ipHD 5 , ipipD 6 , 

etc. These types of quantities can be conveniently simplified by defining an 'electromagnetic radius' D cm ~ 
(ip 2 + H 2 ^J D 6 / M < D. D cnl /D then estimates the fraction of the particle's mass that is of (macroscopic) elec- 
tromagnetic origin. This ratio will come up often, so let 

e ■= < i • (5-6) 

In this notation, (fipD 6 — > e 2 £M, and ipHD 5 — > EM. If <p and H are comparable, the second of these expressions 
is clearly the one that estimates the self-torque. Enforcing (|4.31|) in this case therefore requires 

£ < e 2 . (5.7) 

Conditions weaker than this can be adopted if the particle is dominated by either ip or H a , although it will still be 
true that £ <C 1 (except in special cases). Surprisingly, our approximations place a severe restriction on the charge's 
self-energy. This makes it impossible to ever take the point particle limit in a strict sense (if such a thing is even 
meaningful). Such a procedure is not, however, necessary to answer whether a particle's size is important when it is 
much smaller than the characteristic scales of the surrounding system (the 'physical point particle limit'). 

It is interesting to note that this condition is really independent of any of Dixon's special constructions. In almost 
any slow motion approximation, one would expect that \S\ < \S\T~ 1 and \S\ < MD 2 T~ l . So |*[ ah l| ~ \S\ < e 2 M is 
generic. And this is exactly the condition that led to (|5.7|l . The specific form of the torque needed to show this did use 
Dixon's definition, but the more common one wouldn't have changed anything. This can be verified by substituting 
(|5~T|) into the first term in (|A62|) . 

Also note that when £ ~ 1, it is not clear that any slow motion assumption is even physically reasonable. It 
seems difficult, for example, to avoid the extremely high frequency self-sustaining oscillations discussed by Bohm and 
Weinstein |33f. among others. When £ is small, it is easy to imagine that oscillations like these would have extremely 
small amplitudes, or at least be highly damped in causal systems. But the degree of rigidity required to hold together 
a charge with very large self-energy might prevent this. 



B. Arbitrary Motion 

We now turn back to analyzing the self- forces and self-torques when h a ^ and n a ^ v a . For simplicity, it is 
assumed that the internal dynamics are always 'slow' in the sense of (|5.3[) and 15. 4|) . It is simple enough to relax (|5.5|l . 
so we also allow the first (and only the first) derivative of H a to appear in the self-force. Higher derivatives of both 
cp and H a could be included with relatively little extra effort, although there seems to be little reason for doing so. 



1. Radiative Self-Forces 

At this point, it is useful to separately write out the retarded and radiative fields. Using 1)1. 3[) and H4.17fl - 14.21fl . 

F: c b l{{R) (x)^-^(s)n b \ S ) , (5.8) 

If n a — v a , then this is the same expression found by Dirac from the Lienard-Wiechert potential |12| (we have used 
the opposite sign convention for the field). 

The forces and torques exerted by this field are now very simple to calculate. Inserting (|5.8(l in IJA61J1 . 

* s "iWs) * ~\q 2 n a + o{e 2 £MT- 1 ) . (5.9) 

This is just the Lorentz-Dirac force when n a = v a (see Hl.H ). 

Substituting 1)5. 8|l into l)A62)l . the space-space components of the self-torque are 

J d 3 r pr a /q can be thought of as the separation vector between the 'charge centroid' and the center-of-mass (r a = 0). 
Inverting (|3.31|) , we can clearly convert this component of the self-torque into a 3- vector describing the rate of change 
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of S a . In this form, the (vector) self-torque is just the cross product of the self- force with this separation vector. 
Writing it in this way suggests that it arises due to the self-force not acting through the center-of-mass. 
The time-space components of the self-torque are similar: 

±Jd 3 rH a (r,s)*? R)self rP + o(e 3 £M) . (5.11) 

This clearly vanishes if the 'current centroid' coincides with the center-of-mass. 

Combining all of these results with (|X57|) . ijlH^j) . (|3~33|) . and l|3"31jl . one can find that the particle's motion by 
simultaneously solving 

M ~ -n a 4^ xt + o(e 2 8MT- 1 ) , (5.12) 
S<*P ~ -2^ ] + ^qJd 3 rp(r,s)A a ii® +o(e 3 £M) , (5.13) 

Mn a ~ -* c " xt + ^g 2 n Q + o(e 2 £ MT' 1 ) , (5.14) 
Mv a ~ Mn a + S af3 nf 3 - Mv a - W&hp - 2<i-' a xt 1 e> f) + o(e 2 £ MT' 1 ) . (5.15) 

If the external field varies slowly over E(s) n W, and (SfiH 

can be used to approximate ^" xt and ^^xt • Let 
the minimum characteristic length scale of the external field be denoted by A ^> D, so that \F££ t \ > A|(9F c a x b t | (where 
the absolute value signs are meant to act on each tetrad component of the quantity inside them). It is then convenient 
to assume that 

\>T/£ . (5.16) 

This ensures such that the dipole and higher contributions to the external force are negligible compared to the 
Lorentz-Dirac self-force. 

Without any loss of accuracy, (|5.16|) allows (|5. 1 1(1 - (|5. 15(1 to be written as (making some weak assumptions on the 
magnitudes of the quadrupole and higher moments) 

M ~ -qna&pFgt , (5.17) 
&*> ~ -2e%4 + Q d[a c d d F^) +\qj A 3 r p(r, s)r^n® , (5.18) 

Mh a ~ -qey b F^ t + 2 -q 2 n a , (5.19) 
Mv a ~ Mn a + S a0 np - Mv a - Ae^ b h^ a c F^ t - 2e a a n b & a c F^ t , (5.20) 
where the error terms haven't changed. Similarly, v a can be recovered from (|3.34() : 

Mv a ~ S af} h - 2e a a n b {q\F^ + Q d ^ a c d d F^) -pf d 3 r H (r, s)^r a + o(e 3 £M) . (5.21) 

The assumptions we've made so far - although fairly restrictive - are clearly not sufficient to recover the standard 
point particle result. Still, there must exist some class of charges which do behave in this way, and it is interesting to 
characterize it. Before doing so, it must first be mentioned that the s in l|l.l|l is a proper time, while s is not. It is 
close, though. Letting s = s(s) and s'(s) := ds/ds, 

' (5.22) 



\fv a v a 

~ l + o(e 4 ) . (5.23) 

The fractional difference between the center-of-mass 4- velocity and v a is therefore of order e 4 . In particular, the 
triad components of this 4- velocity will differ from v a by terms no larger than e 6 . Given the error estimate in 1(5. 21j) . 
these differences are negligible whenever £ > e 3 . 

The proper acceleration can now be written in terms of v a and v a : 

z" a = s"v a + (s'fv a . (5.24) 
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It follows from 1(5.20(1 and (|5.22() that s" ~ e 4 /T, so z" a differs from v a by terms of this same order. (|5.15|l then 
implies that these two quantities are interchangeable whenever £ > e 2 . 

Applying similar arguments, our results can be easily compared to the Lorentz-Dirac equation by writing it in the 
form 

(Mv a ) hD ~ -gey b F^ t + ^q 2 n a + of^MT' 1 max(e,£)) , (5.25) 
(Mn a v a ) LD ~ -qn a e^v F^ t + ofe^MT^maxie^)) . (5.26) 

In order to avoid overcomplicating the discussion, we shall say that our equations of motion reduce to the Lorentz- 
Dirac result if Mv a and Mn a v a match (15.25(1 and ((5.26(1 up to terms of order e 2 MT~ 1 max(e 2 , £). Noting that 
n a v a = —h a v a , (|5.1 9(1 shows that the the 'temporal component' of the body's acceleration always matches 1(5.2611 to 
the required accuracy. 

The same is not true for the spatial acceleration. This shouldn't be too surprising, though, as the Lorentz- 
Dirac equation was never intended to describe spinning particles. S ab should therefore be set to zero (at least 
instantaneously) before any reasonable comparison can be made. Even this isn't quite sufficient, though. The 
situation can be remedied by assuming that 

\S\ < eM J Dmax(e 2 ,£) , (5.27) 
Q^ifeef < e 2 Mmax(e 2 ,f) , (5.28) 



along with a similar restriction on 



Q[a pb]c 



oxt e a n b 



. Then 



\v\ < e 2 max(e 2 ,£) , (5.29) 
\M\ < e 2 MT" 1 max(e 2 ,£) , (5.30) 
\S\ < e 2 Mmax(e 2 ,£) . (5.31) 

If £ > e 2 , M vanishes up to the maximum order that we can calculate it. The same is not necessarily true for S a ° and 
v a , although they remain sufficiently small that ((5.25|) can now be recovered to the desired accuracy. There therefore 
exists a regime in which the equations of motion derived here reduce to the Lorentz-Dirac equation. Considering z a 
to be the only observable (as would be reasonable for an extremely small particle), this completely recovers the usual 
point particle result. 

While this conclusion is not particularly surprising, it is interesting to note how restrictive the required assumptions 
are. ((5.28(1 is particularly difficult to satisfy. For example, when the force is approximately given by the Lorentz 

(monopole) expression, external field magnitudes up to ~ M(|g|T) 1 are allowed. In fields this large, 15.28(1 implies 
that the magnitudes of the dipole moment must be less than ~ e£qD (when £ > e 2 ). This is an extremely limiting 
condition in cases where the self-energy is small. Given (13.15(1 . there do not appear to be any rigid (ip = H a = 0) or 
nearly rigid charge distributions that could satisfy it. One either needs to choose a very special class of charges, or 
considerably restrict the maximum allowable field strength. 

Before continuing, it should first be mentioned that our definitions of the charge and three-current densities are 
slightly unusual. These quantities would usually be defined with respect to an orthonormal tetrad adapted to z' a (s(s)) . 
We have instead defined these quantities in terms of a tetrad with temporal component n a . Translating from one of 
these definitions to the other would only introduce fractional changes of order e 2 , which are usually irrelevant at our 
level of approximation. 



2. Retarded Self-Forces 

It is generally accepted that detectors placed outside of W will measure F^ t + F^ { ,s rather than F^ t + F^ { , R y 
But infinitesimal elements of an extended body cannot 'know' that they are part of a larger whole, so they must 
couple to this same field. It is therefore reasonable to consider only the retarded self-field 'physical.' Recalling that 
this is equal to the sum of the radiative and singular self-fields, any relevance of F^ l{ ^ itself should be derived by 
showing the effects of the singular self-field are irrelevant in certain cases. 

To this end, we now examine the (presumably) realistic case in which the particle interacts with its full retarded 
self-field. Essentially all of the steps in the previous section carry over identically in this case, although most involve 
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considerably more calculation. Starting with (|3.7(l . 13.8(1 . and (|4.17|) - (|4.21|) . the retarded field can be shown to be 
(dropping the '— ' subscript) 

F&(x) * 2 J dV ± L [-(r - r>Ye [ > b] (l + \n f \r r') - ±(^(r - r') f + |i? 2 |^| 2 

+ -h p Wrp(r - r') 7 ) + -R 2 {r - r') a e [ £^ - -R 2 n^ a n^ 1 1 + ^(Sr - r')A 

2 / 2 2 v 2 y 



+ (r"V) ) (r - r') Q eL a e^ (\ + in^r + r% + ^(n 7 (r + r')- 



+ ii? 2 |n| 2 - -jVT?VV 7 r^ + n^ejgii 2 (-n<{r + r') 7 + -i?H 2 - R- 2 v~<(r - r'\ 



+ ip{r - r') a e^e b yP - #V a eJ ] .R 2 j . 

The self- force is now obtained by inserting this into (|A61fl . which results in an expression of the form 

KM - - Jd 3 rd 3 r'F a (r,r\s) , 

r 1 (T a {r,r',s)+T a {r\r : s)} , 



d 3 r dV 



(5.32) 

(5.33) 
(5.34) 



where T a {r 1 r' , s) represents the force density exerted by a charge element at r on a charge element at r'. If Newton's 
third law were correct, J :a {r,i J ,s) = —J- a (r',r,s) (so ^^if would vanish). Of course, this is does not quite hold for 
the electromagnetic field (or any other fully observable field), so there is a nonzero self-force. ((5.34(1 therefore gives a 
precise form to the intuitive idea that self-forces measure the degree of failure of Newton's third law. 

After removing the components of J- a which reverse sign under interchange of r and r' , it can be shown that 



Mn° 



¥£t + -q 2 n a - M cm h a - J d 3 rd 3 r' — jp(r, s)H f3 (r', s)R 
- iT"V(r - r')^j - (r - r') a h + p(r, s)H f3 (r' , s)R 2 + v^dyfr^r, s)) H a (r' ', s){r - r'f 



tp(r, s)H^(r ,s)(r — r ) Q r 7 — ( (p(r, s) 



An\r\ 



i# 7 (r', s) (R 2 n a r 7 - R 2 S"h 



r 



+ ((r - r')% - <5«(r - r'frp)n"{r - r% j + H J (r' \ s) ((r - r') a r 7 - 5«(r - r')%) 

+ o(e 2 £MT- 1 ) , 
where it was useful to define an 'electromagnetic mass' 



(5.35) 



d 3 rdV 



+ <p(r,s) 



- - 
i? 3 | 2 

9 



p(r, s)p(r', s) - ifyfa s)H /3 {r', s)(l + n> '(r + r')- 



47r|r| 



p(r',s) 



(r - r')% f 1 + in?(r - r') 7 ) - Ji? 2 n% 



(5.36) 



The form of ((5.35(1 suggests an effective inertial mass m := M + M em . Although M will very rarely remain constant, 
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m is often conserved, or at least varies slowly. To see this, substitute (|5.32|l into (|A61fl and use (|3.32() to show that 



d 3 rdVi)(p( M )- " 



An \ r • 



p(r', s)(r - r')% - ~p(r', s) (i? 2 ^ 



-n Q (r-r%(r-r')%) -#V, S ) Qi? 2 (n>% + r>| 2 ) - (r - r') Q (n Q r. 
- r a n,) (l + in ff (r + r%)) - H fJ (r)H^r')R 2 - l -H rj {r)H^ \r')R 2 rO '(r + r')- 



p{r',s) 



R 2 H (r,s)h ^l + in Q (3r-r')a^) - « 7 d 7 (w(r, s)) (r - l+o(e 2 £MT- 1 ) . (5.37) 



Combining this with l|5.36[) and the definition of m, 

if/,. 7 



d 3 rdV 



R 3 



<f(r,s) - 



Att \r • 



F 7 (r',s) 



(n 7 (r - r')°V Q - r 7 h a (r - r')a) 



(l + ~n A (r + r') A ) + ^i? 2 (N% + n%n 7 ) 
i; 7 d 7 ( W (r, s)) (r - r')' 3 | + o(e 2 £Af/T) , 



•p(r',s) 



R' z H /3 (r,s)n p ( 1 + in a (3 



(5.38) 



which is considerably simpler than (|^ 

Although we now have all of the results necessary to compute each tetrad component of W^j! to second order, 
the resulting expressions are extremely lengthy (and correspondingly difficult to interpret). It was also seen when 
examining the radiative self-fields that the second-order terms in the self-torque were not necessary to reach the point 
particle limit. For both of these reasons, we only compute it here to first order. Combining (|A62|) and l|5.32[l . 



a/3] 

self 



/ 



d 3 rdV ■ 



,[«t-/3] 



R 



o P( r i s )p( r ', s ) + H i( r > s)i? 7 (r', s) + <p(r, s) 



A-rtM 



p(r',s) 



+ o(e 2 £M) , 



(5.39) 



n <* e b *sclf 



i Jd 3 rd 3 r'^ ^p(r,s)H a (r>,sye (r - r') a (l - \w{r - r') 7 ) + l -R 2 n 

- (<p(r, s) - -fr] R 2 r 7 h^H^(r',s)+2r 7 (r-r') V H^(r',s)(l + ln^(r + r') 
\ Am \t\ ) y V 2 

- r a H a {r f ,a){r p n'r' 1 - r //3 n 7 r 7 ) + (r - r') a r a {h p {t\ s)n 7 r 7 - r f3 h^H' f (r', sfj 
+ o(e 2 £M) . 



(5.40) 



VI. SPECIAL CASES 



A. Current-Dominated Particles 



The full equations of motion derived above are obviously quite complicated,, so it is useful to specialize the discussion 
somewhat. First assume that the particle's charge density is sufficiently small that it can be entirely dropped without 
losing any accuracy. This requires that ip < e 2 H. Then (|5.39(l becomes equivalent to 

S<*P ~ -2*["f + / d 3 r dV — — Ky(r, s)H^(r', s) + o(e 2 £M) . (6.1) 
This can be quite large. In order to make sure that it satisfies 1)4.32(1 . let 



£ < e 



(6.2) 
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Now, M cm < eM, so the remaining equations of motion can be derived from Q5.350 . (|5.38(l . and (|5.40() : 

m ~ -« Q *e X t + o(e 3 mT- 1 ) , (6.3) 

mn a ~ -*^ xt + o(e 3 mT- 1 ) , (6.4) 

mv a ~ S ,Q/3 n /3 - 2*J, < $e> 6 + o(e 3 m) . (6.5) 
An expression for v a is obtained by differentiating (|3.34() . which gives 

d 3 r dV ^ ^ g 7 (r, s)iF(r', s) + o(e 3 mT- x ) . (6.6) 

The presence of a (potentially) significant self-torque here is completely different than the situation that arose when 
considering only the radiative component of the self- field. Even adopting conditions l|5.16l) . I|5.27l) . and l|5.28[l would 
not generically recover the Lorentz-Dirac equation. 

Instead, we find (assuming £ > e 2 ) that (|5.29|l and (|5.3U|I (with M — > m) would remain correct, although (|5.31|l is 
weakened to IS"! < e 2 m. In this case, the center-of-mass line is governed by 

mv a ~ -qey b F^ t + -q 2 fi a + 2M em D l ^v^v + o(e 3 mT- 1 ) , (6.7) 

where 

ms) := m~ / dV d3r ' (5) H ^ s)jff7(r '' s) • (6 - 8) 

Given Ij5.36|l . _Dg appears to be related to the shift in the position of the effective center-of-mass due to the electro- 
magnetic self-energy. 

The force that it generates clearly becomes irrelevant when M em D^v^V/3 < e 3 m/T, which occurs when £)g (nearly) 
coincides with i) a , or more generically if 

\D H \ < eD . (6.9) 

This happens, for example, in cases where H J (r, s) ~ ±iJ 7 (— r, s). I|6.9|l can therefore be thought of as a restriction on 
the allowed asymmetry in the particle's current structure about the center-of-mass. Given Ij3.15|l . a similar intuitive 
interpretation can also be applied to (|5.28JI . so there is some overlap between these two conditions. It does not appear 
that either one strictly implies the other, however. 

Also note that if (|6.9|l holds, (|6.2() is no longer required. The case £ ~ 1 is then quite interesting if the spin and/or 
dipole moment are non-negligible. When this happens, the left-hand side of (|6.5|l needs to be replaced by Mv a , which 
means that M and m must both be kept track of. The expression for M is not simple, so this is a considerable 
complication. 



B. Charge-Dominated Particles 



Treating the opposite case, we now consider particles with very small internal currents. In particular, let H < e 2 ip. 
Given (|A40|I . this is about as small as H could possibly be without fine-tuning. I|3.8I) and 1)5. 3|) show that the internal 
currents which arise due to charges rearranging themselves (via elasticity, rotation, etc.) are bounded by this same 
amount. So we actually have that \j\ < e 2 ip. 

(|5.39|l now reduces to 



-2* 



d r d r — — — p(r , s) 
R 



p(r, s) + 2 ip(r, s) 



Am \r r 



o(e £m) 



Once again, this expression can easily violate l|4.32|l unless £ < e. Assuming this, 



m ~ -n a *° xt + I d 3 rd 3 r' -v a d a p(r,s)p(r',s) +o(e 3 mT- 1 ) , 



cxt ,- -q<n a + o (e*mT- 



mv a ~ 5 Q/3 n /3 -2^ x b t 1 e> 6 -|-o(e 3 m) 



(6.10) 

(6.11) 

(6.12) 
(6.13) 
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By analogy to the current-dominated case, we can now define 



D a ■= — 



d 3 rdV(-)p(r', S ) 



p(r, s) + 2 1 ip(r, s) 



<1 



47r|?'| 3 



(6.14) 



Again, this looks like the center-of-electromagnetic mass, although the interpretation isn't quite so direct as it was 



for L>g. 



Now let £ > e 2 , and adopt l|5.16[) . I|5.27|l . and Ij5.28|l . As before, these conditions imply that \v\ < e 3 , so m < e 3 mT 1 . 

The Lorentz-Dirac equation is recovered when M em D^v^vp < e 3 mT _1 . If the acceleration is not restricted to be in 
a specific direction, then it is convenient to satisfy this condition by requiring 



\D V \ < eD 



(6.15) 



This clearly holds when (for example) (p(r, s) ~ ±ip(—r,s), so \D V \/D can be considered a measure of the charge's 
internal symmetry. 

The two examples discussed here are by the far the simplest, although other ones may also be of interest. For 
example, the condition \j\ < ep could have replaced \j\ < e 2 p in this section. This would allow for some charge- 
current coupling without introducing undue complexity. Similarly, one might be interested in the case where p> < eH , 
or the general scenario where p> and H a are symmetric (or antisymmetric) about r = 0. 



VII. CONCLUSIONS 



A. Results 



Starting from the fundamental equations (|2.1|) - (|2.4|) . we have now derived equations of motion describing a very 
wide variety of classical extended charges in flat spacetime. One of our main motivations has been to investigate the 
validity of the commonly-held notion that 'small' charges can often be treated as though they were perfectly pointlike. 
Because of this, the precise implications of each approximation have been emphasized, and all of our results have been 
kept as general as (reasonably) possible. 

To review, charge distributions were considered where all significant length scales remained of order D(s) within 
each time slice S(s). An acceleration timescale T 3> D was then defined to place lower bounds on \fi\~ 1 and |ri|~ 1//2 . 
This was assumed to relate to the timescales of the body's internal motions as well. Specifically, the s-derivatives of the 
'charge potential' p> were required to satisfy (|B . 3|) and (|5.4(l . Similar restrictions were placed on the derivatives of H a 
as well. The last major restriction required that the spin and torque play a relatively small but non-negligible role in 
the system's behavior. This was formalized by assuming that \n — v\ = \v\ remained less than about e 2 = (D/T) 2 <C 1. 
All of these assumptions are kinematic, and cannot automatically be assumed to hold without any regard for the sizes 
of the external fields, spin, and so on. Self-consistency is preserved by requiring the angular momentum, force, and 
torque obey the bounds discussed in Sec. IIV Al 

The restriction given on S ab essentially states that the body's rotational period can't be less than T. The relations 
satisfied by the force and torque are not quite as simple, though. Both of these quantities depend on J a , so any bounds 
placed on them must affect the class of allowable current distributions. While this can be taken into account in different 
ways, we have chosen to impose restrictions on the dipole moment and the 'relative self-energy' D cm /D — £. If the 
dominant force is just —qF^Vb, Q ah must satisfy (|4.37() . 

The bound on £ is more complicated to summarize, as it is sensitive to the specific case under consideration. 
The calculation using the radiative self-field is the least restrictive; allowing all £ < 1. One could technically let 
the electromagnetic radius exceed the physical one in this case, although we take the point of view that this would 
be too unphysical. Regardless, taking into account the full retarded self-field shows that there exists a significant 
class of charge distributions where the self-torque will become too large whenever £ ^> e 2 . In scenarios where the 
charge-current coupling is negligible, the relative self-energy can be of order e. Still, there are certain special cases 
such as spherical symmetry where the restriction on £ can be completely relaxed. 

Despite appearances, these assumptions are really no different than the standard slow-motion approximation. 
Using it, the radiative and retarded self-fields in a neighborhood of W were shown to be approximated by l|5.8|l and 
(|5.32l) respectively. These expressions are quite general, and although the retarded field is written using the peculiar 
constructions of Dixon's theory, it can easily be translated into a more conventional notation by using 7|) and (|3.8|) . 
(|4.17l) - (|4.21() also provide a convenient starting point for this. 

These expressions for the fields were then combined with the exact forms of Dixon's equations of motion to yield 
the results contained in Sec. All of the approximations adopted here were therefore applied only to compute the 
fields. In the (unrealistic) case that the body was assumed to couple only to the radiative component of its self-field, 
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the center-of-mass line was found to evolve according to H5.12[) - <|5.15[1 . The external forces and torques that appear 
in these expressions are given exactly by l|A61|l and (|A62I) . or more intuitively by (|A63(I and l|A64|l . 

In order to investigate the Lorentz-Dirac limit, it was first necessary to require the particle's radius to be sufficiently 
small that the dipole and higher contributions to the external force could be considered negligible without having to 
throw away the Lorentz-Dirac component of the self-force. This idea was summarized in the condition (|5.16l) . which 
transformed the equations of motion in the radiative case into (|5.17(l - (|5.21(l . Recovering the Lorentz-Dirac equation 
was then seen to require bounding the angular momentum and dipole moment by (|5.27|l and l|5.28|l respectively. 

Since this limit is meant to describe a non-spinning particle, the bound on the angular momentum is hardly 
surprising. (|5.28|l is rather different, though. It shows that even when ignoring the 'singular' portion of the self-field, 
the Lorentz-Dirac equation does not apply to all small charge distributions. If a problem were to involve very large 
(external) field strengths and small self-energies, the class of charges which move like a point particle would in fact 
become extremely small. Oddly enough, this problem disappears for large self-energies. In these cases, the restriction 
on the dipole moment is no stronger than was required for self-consistency of the initial slow-motion assumptions. 

This issue is made considerably more complex when the body is allowed to interact with its full retarded self- 
field. Without placing any restrictions on the spin or dipole moment, n a was found to be given by (|5.35|1 . This 
involves an effective mass m = M + M em , which involves an electromagnetic contribution given by l|5.36[l . It does 
not generally remain constant, but rather evolves according to l|5.38[l . There are also self-torques given by (|5.39|) 
and (|5.4U|) . Combining these expressions with (|A57I) and 13. 3411 recovers the full equations of motion. These are our 
central result. They provide a considerable generalization of the Lorentz-Dirac equation. 

They are also much more complicated than the equations derived using only the radiative portion of the self-field. 
The majority of this complexity arises from interactions between <p and H a , which - given l|3.7|) and (|3.8|) - can 
be viewed as charge-current couplings in the center-of-mass frame. Predictably then, our results are considerably 
simplified when only one of these quantities is significant, ip will (generically) drop out of the equations of motion 
only if it is of order e 2 H or less. The reverse is also true. Both of these cases were discussed in detail in Sec. IVII 
It should be emphasized that these restrictions on the relative magnitudes of the charge and current densities are 
exceptionally strong. Simply saying that one of these quantities is much larger than the other without any further 
qualification is not sufficient to remove the coupling terms. 

Even in these cases, though, the previous restrictions on A, S ab , and Q ab must be supplemented by either l|6.9|) 
or I|tj.l5fl m order to recover the Lorentz-Dirac equation. Intuitively, these relations ensure that the 'center-of- 
electromagnetic mass' is not too far from r = 0. This is roughly what was already required by (|5.28|l . although the 
two conditions are not quite the same. 

Note that the assumptions we have given to be able to derive the Lorentz-Dirac equation are not exhaustive. They 
were the most obvious choices obtained by examining our general equations of motion, but are slightly more restrictive 
than necessary. Still, they seem to be reasonably effective for most systems that are not finely-tuned. 

B. Discussion 

In summary, it was shown that the Lorentz-Dirac limit is rather delicate, and that the radiative self-field is rarely 
an adequate replacement for the retarded one. It is now interesting to speculate how these results might generalize 
in curved spacetime. One might be interested in the motion of a charged particle moving in a background spacetimc, 
an uncharged body allowed to generate its own gravitational field, or even the general case of a massive charge. With 
the exception of this last possibility (which doesn't necessarily follow from the two simpler calculations), such systems 
have been considered in the past by a number of authors |4l Ha . fl6t 1171 l25l 12a . l36| . All of these calculations have either 
applied a point particle ansatz, or considered only very special (though usually astrophysically motivated) classes of 
extended bodies. Both such methods have agreed with each other, although it is clear that some extended bodies 
must exist which do not move like point particles. 

Finding the scope of the existing MiSaTaQuWa equation 01 ~ as we ll as its generalization - would likely require 
an analysis similar to the one given here. The methods used in this paper have in fact been chosen specifically for 
their ability to be easily applied to fully dynamic spacetimes (except for those used to expand the self-field). To see 
this, it should first be mentioned that all of the advantages of Dixon's formalism are retained without approximation 
in full general relativity |23l |24| . 

This provides two essential results. The first of these is a natural notion of a center-of-mass line. This definition 
has been proven to satisfy nearly all of the intuitive notions that one might expect of such an object |29j . so it can be 
considered a reasonable measure of a particle's 'av erag e' position. Just as importantly, it is exactly determined by a 
finite number of ordinary differential equations p^.l28j|. The analogous equations here were l|3.32[) - (|3.34[) . and l|A57[) . 
These expressions are only slightly complicated by the transition to curved spacetime. 

The second contribution of Dixon's formalism is the set of stress-energy moments itself. These are all contained in 
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the 'stress-energy skeleton' T ab (r, s). Here r is a coordinate in the tangent space of z(s), which is an essential point 
that was not obvious in flat spacetime. Importantly, the constraints on T ab implied by the generalized form of i|2.4|) 
are exactly the same in all spacetimes. This obviously includes the fiat case summarized in the appendix. Just as 
(|A24|) gave the general form for any current skeleton J a satisfying the constraints, one can also find all possible forms 

of f ab 113. 

As with current skeleton, not all of these possibilities are physically reasonable. Some will be singular, and others 
will have supports extending to spatial infinity. We conjecture that such cases can be removed exactly as they 
were in the appendix for the current moments. This would leave a set of 'reduced moment potentials' that could 
be arbitrarily generated from some simple recipe. An automatic byproduct of this reduction process would be a 
relatively straightforward method of generating T ab from the potentials. The analogous results in this paper were 
that the current skeleton was determined by the freely-specifiable functions cp and H a ^ . These in turn generated J a 
via ipT7|l . (23, and (|3~11fl) . 

In the end, such reduced potentials for the stress-energy moments would serve (in combination with solutions of 
the ODE's for the linear and angular momenta) as a nearly background-independent way of specifying physically 
reasonable conserved stress-energy tensors. Besides the intrinsic elegance of such a construction, it would also solve 
several problems that did not arise in the present paper. The most important of these is the location of the center-of- 
mass. At any s, this is clearly determined by T ab on S(s). The definition is highly implicit, however, and it is almost 
always impractical to apply in practice. Specifying the matter via T ab would avoid this problem. Constructing it in 
the natural way would start with the center-of-mass position. It would also incorporate p a and S ab automatically. 
These quantities obviously must solve certain differential equations, and writing down a T ab on each time slice with 
the correct momenta would be extremely difficult with any generality. These issues were the main motivations for 
discussing the current moments in so much detail in the appendix, as well as the use of ip and H a throughout this 
paper instead of J a . These choices slightly obscure the results here, but would be essential in any generalization. 

With this formalism in place, one would then need to compute the metric to find the body's motion. As in the 
electromagnetic case, this is the most difficult step. Of course, Einstein's equation is considerably more complicated 
than Maxwell's, so it may be prohibitive to carry out the calculation by hand (except in special cases). Regardless, the 
above procedure could be integrated into existing numerical relativity codes. This would then allow one to rigorously 
study spacetimes with nonsingular matter fields without having to solve the conventional elasticity (or Navier-Stokes) 
equations. The most convenient such systems in this formalism may not represent the most astrophysically interesting 
types of matter, although there is no shortage of important problems in numerical relativity where the details of the 
matter distributions are not a primary concern. 

Whatever the results of such inquiries, our results in electromagnetism can be used to speculate how point particle 
methods might break down in gravitational self-force problems. For this, it is useful to give an intuitive explanation 
for such failures in Maxwell's theory. The simplest of these derives from the fact that the center-of-mass does not 
generally correspond to anything that could be called a 'center-of-charge.' In these cases, electromagnetic self- forces 
effectively act through a lever arm. This induces a torque, which in turn affects the particle's overall motion. The 
situation is made particularly complicated when the self- force is broken up into several pieces (as is natural). Each 
such piece tends to act through a different point, and there's no particular reason that any of them should coincide 
with the center-of-mass (although the initial restriction (|4.37|) on the magnitude of the dipole moment does help) . 

In contrast, one might expect that gravitational self-forces would always act through the center-of-mass; effectively 
removing this type of effect. This isn't necessarily true, though. The center-of-mass definition only requires that 
t abc n n c = 0, which does not necessarily mean that t abc completely vanishes. The portion of the self-force acting 
as an effective mass could also produce a significant torque. (|A47|) shows, however, that t abc does vanish whenever 
gab _^ q This is the only case in which the MiSaTaQuWa equation can be reasonably expected to hold, so imposing 
it at the start should remove any problems. Since the dipole moment here is entirely dependent on S ab (unlike in the 
electromagnetic case) it should be relatively easy to arrange for an initially-vanishing angular momentum to remain 
small. These sorts of effects might therefore be expected to be less troublesome in gravity than in electromagnetism, 
although they will probably still exist. 

This type of mechanism did not lead to most of the complications found in this paper, however. These were instead 
due to couplings between the charge and 3-current densities as viewed in the center-of-mass frame. In gravity, the 
situation could be even worse. Similar interactions might exist between mass, 3-momentum, and stress densities. 
Intuitively, though, it would appear that adopting appropriate energy conditions should considerably soften these 
interactions. 
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APPENDIX A: DIXON'S FORMALISM 
1. Current Multipoles 

This portion of the appendix reviews Dixon's decomposition of the electromagnetic current vector into multipole 
moments. It also derives natural 'potential functions' that can be used to generate sets of moments for all physically 
interesting current vectors satisfying (|2.3() . These objects are shown to determine J a in a simple way. The main 
reasons for these constructions are explained in Sec. IV11 hSI although some secondary points are also mentioned in 
Sec. Mil In short, analogs of these steps would become essential in the gravitational self-force problem, so we include 
them here to allow a relatively straightforward generalization. 

The notation here will be that defined in Sec. Mil Using it, we can define multipole moments for the current vector. 
Such objects are usually constructed by integrating a source function against a suitable number of radius vectors. 
Defining r a = x a — z°(s) for x 6 S(s), one might therefore expect that for n > 0, the 2n-pole moment could be given 
by 

Q«-i-*»«( a ) : = Jd 4 xr bl ■■■r K J a (x- z(s),s) , (Al) 

where we have assumed the existence of some distribution J a (x 1 s) which is related to J a (x), but has compact support 
in x. This will be called the current skeleton. Note that (|A1|) automatically implies that 

Q bl - b " a = g( fc i- b ") a (A2) 

for all n > 1. 

It is rather cumbersome to keep track of each Q " directly, so we instead define a generating function 

G a (k, s) := Q a (s) + V y -^-k hl ■ ■ ■ k bn Q b i- Ka (s) . (A3) 

n—l 

An arbitrary multipole moment can now be extracted from this in the usual way: 

Q^- Ka (s)=i n (d bl ■■■d b "G a (k,s)) . (A4) 

V / k=0 

G a is therefore completely equivalent to the set {Q a ,Q ba , ■ ■ ■}■ 

Although this is a useful property, the definition of G a is not simply a mathematical convenience. Using (|A1|) and 
(IA3I) , it can be shown that 

G a {k lS ) = Jd 4 rj a {r,s)e- tkr . (A5) 

The Fourier transform of G a , 

G a (r,s) := [d 4 kG a (k,s)e lk r , (A6) 



is therefore proportional to J a : 



G a (r,s) = (2ir) 4 j a (r,s) . (A7) 



This shows that if J a and J a equivalent in an appropriate sense, the set of moments can be used to completely 
reconstruct the current vector. 

In order to relate the current to its skeleton, it is convenient to think of J a as a linear functional on the space of 
all C°° test functions with compact support (as is typically done in distribution theory). In particular, knowing 

(j a {r,s),<l> a {x)) := Jd 4 xJ a (r,s)Mx) (A8) 

for all suitable test functions <fi a (x) can be used to define J a . An analogous statement can also be made for J a . These 
two objects can then be related to each other by writing (J°, <f> a ) in terms of ( J a , </> a ). The latter expression depends 
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on s, while the former one does not. It is therefore most straightforward to link the two by simply integrating out 
the s-dependence: 

(j Q (x),0 a (x)) = Jds (j a (r, S ),Mx)) , (A9) 

= Wf! dS ( GQ (M,0a(fc)e-^ (s) ) . (A10) 

Following [22^. we take this (along with (|A3f> ^ to define what is meant by saying that that the Q "'s are 'multipole 
moments of J a .' 

This is not uniq ue definition, however. To remove the remaining freedom in a useful way, we simply state Dixon's 
results [iSlllllJ. Let 

n bl Q bl - b - l[b " a] = 0, (All) 
Q^- h ^ = , (A12) 

for all n > 2. Also assume that the monopole moment has the special form 

Q a = qv a , (A13) 

where q is the total charge as it is usually defined. 

Now choose a test function of the form <j> a (x) = d a <j>(x), with <j)(x) itself also a test function. Then l|A4(l . (| A10|) . 
and l|A12fl can be used to show that 

<a Q j a » = -(J^dacjy) , 

= 7A4 J&s(k a G a {k, S ),mz- ik - z[s) ) , 

= - Jdsqv a ( S )d a ^(z( S )) , 

dsq±<j>(z(s)) . (A14) 

This must vanish for all 4>, which can be ensured by simply requiring that 

= 0. (A15) 

This (trivial) evolution equation is the only one implied by 1)2.3)1 . It can be shown that moments satisfying l)A2)l . 
(|A11|) - (IA13|I . and (|A15|) describe any J a with the given prop erties in a uniquely simple way [20I I22I l23j. A precise 
statement of the theorem that was proven is contained in [23j. 

These conditions allow a great deal of freedom in choosing different moments. The fact that there is only one 
evolution equation implied by 1)2. 3f) does not mean that the higher moments necessarily remain constant. Rather, 
they can be given an almost arbitrary time dependence. This is a reflection of the fact that we have not yet chosen 
to model any particular type of matter. Specifying how the higher moments change in time is essentially equivalent 
to choosing an equation of state. Although this must be given in order to have a well-defined initial value problem, 
self-forces and self-torques can be written down without ever having to explicitly evaluate the time derivatives of the 
current moments. This allows us to derive equations of motion valid for a very large class of systems. 

In order to do so, we first need to pick out sets of moments (or equivalently their generating functions) which 
represent physically reasonable current vectors. Despite appearances, this cannot be done arbitrarily. Dixon's theorem 
ensures that any nonsingular current vector with support W can be described by a set of moments with the given 
properties, although the reverse is not necessarily true. Extra conditions need to be imposed in order to ensure that 
the J a associated with any particular G a (or J a ) has the correct smoothness and support properties. 

To gain some insight into this, fix some test function <j) a . From this, construct a second test function (j)' a which agrees 
with (j) a everywhere except in an infinitesimal neighborhood of Z. Assume that the support of <f>' a does not include Z 
itself. (| A3|) and i|A10|) then show that any finite number (and only finite number) of multipoles can be changed without 
affecting (J a ,(j)' a ). The same cannot be said for (J a , O ). But for any physical current vector, (J a ,<f)' a ) ~ {J a ,cj) a }. 
This shows that any finite subset of an admissible collection of moments is completely determined by its complement. 

Because of this, it is not reasonable to impose conditions directly on the individual moments to ensure that their 
associated current vector is physically acceptable. Such restrictions are most easily stated in terms of G a or J a , but 
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doing so first requires finding how J a is affected by (p^TTll - (fXT3|) . and This is now done by finding the general 

form of G a , and then taking its Fourier transform. 

It is shown in pflj that the constraint equations are precisely equivalent to requiring that the generating function 
have the form 

G a (k,s) = G a {1) (h-k,s) + (n- k)G a {2) (h-k,s) , (A16) 

where n ■ k := n a k a , (h ■ k) a — h b k b , and n a G? 2 <. = 0. Also, 

fc Q G a (fc,.s) = qk a v a (s) , (A17) 
d [a G^ 2) {h-k,s) = . (A18) 

It is now useful to take Fourier transforms of these equations to find their equivalent forms when representing the 
moments by J a (= G a /(2ir) 4 ). doesn't take on any special form, although (|A18|I shows that 

(r c / 1 [ Q G^(r, S ),0(r))-O. (A19) 

This equation is solved by any G^ of the form 

Gf 2) (r,s)^h a b (s)r b G {2) (r, S ) , (A20) 

for all functions G(2) (Despite the notation, it will shortly be clear that the inverse Fourier transform of G(o\ does 
not exist in general.). Note that this is not the most general solution of l|A19(l . A term of the form g a (n ■ r, s)S 3 (h ■ r) 
may also be added to G%y although we choose to ignore this possibility. 
In any case, (|A16|) and l!A20|) give an explicit form for G a 

G a (k,s) = G a {1) - ih a b r b n c d c G (2) . (A21) 

Since G"^ and G° 2 ^ are independent of n ■ k, their (four-dimensional) Fourier transforms must be proportional to 
5(n ■ r). It is therefore natural to define quantities A, B a , and G such that 

G a {1) (r,s) = (2vr) 4 (5(n • r) (A{h ■ r, s)n a + B a {h-r,s)^ , (A22) 
G( 2 ) (r, s) = -i(27r) 4 (5(n • r)C(h ■ r, s) , (A23) 

where B a n a — 0. Combining these expressions with (|A21|) . J a is found to have the form 

J a (r, s) = S(n ■ r) (aQi ■ r, s)n a (s) + B a (h ■ r, sfj - S'(n ■ r)C(h ■ r, s)h1,(s)r b . (A24) 

This is further restricted by l)A17(l , the Fourier transform of which becomes 

(daPfas),^) = (2ir)- 4 i(k a G a (k,s),0(k)) , 

= -(27r)-V a (s)(T,9 a 0(r)) , 

= ^°(s)(d a 5 4 (r),<Kr)) • (A25) 

Comparing this to the divergence of (|A24(I shows that 

q5 3 {h-r) = A(h-r 1 s)~h a b d a (r b C{h-r,s) S j , (A26) 
qv a d a 6 3 (h ■ r) = d a B a (h ■ r, s) . (A27) 
The solution to the 'homogeneous' analog of (|A26|I . 

h a b d a (r b C (a) (h-r,s)) =-q6 3 (h-r) , (A28) 
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C(o){h-x,s) = - 



47r|r| 3 



where |r| 2 := —h a i,r a r b > 0. For future convenience, we now define new functions tp(r,s) and qo(s) such that, 



C = N [ ip — 



ffo(s) 
47r I r 1 3 



(A29) 



(A30) 



N is just the lapse, as given in l|3.5fl . 

In terms of qo and ip, l|A26|l now becomes 



A(h-r,s) = h a b d a [r b Np(h-r,s)) + (q - g ( s )) S 3 (h ■ r) + 



go TTTq 

47r Irl 3 



(A31) 



J a is therefore given by l|A24|) with A and C having the respective forms (|A31|) and i|A30() . We also have that B a 
satisfies (JA27|1 . This effectively parameterizes all sets of moments satisfying Dixon's constraints. 

I|A9J) can now be used to relate J a to J a . It is convenient to do this in the (r a , t) coordinates defined in Sec. 1 1 1 1 1 
Using them, l|A9(l becomes 



(A32) 



(J a ,<M = J dsj d 3 rl dTN(x)J a (x- z(s),s)Mx) ■ 
Applying (|A24ll . the r- integral in this equation can be carried out explicitly: 



ds / d A r 



An a + B a + N- 2 e a a r a h,3V f3 C + N^e^v 13 dp^C) 



N- l Ce a a r«^ 
as 



(A33) 



where all quantities are evaluated at s. 

Commuting the s-integral with the spatial ones and integrating the last term by parts, 



{J a Aa) = d 



ds 



A + -hpr? 



B a + v^dg r a 



C_ 
N 



d (C 



ds \N 



(A34) 



But we also have that ( J a , (f> a ) — J d 3 r J ds NJ a cj> a , which gives us an obvious way to explicitly write J a in terms of 
A, B a , and C. 

The charge density with respect to the tetrad frame takes on the particularly simple form 



p := n a J a , 

= d a (r a ip) +N- 1 (q-q )5 3 (r) 



(A35) 
(A36) 



Physically, p(r a , s) must be nonsingular and have support W. Without loss of generality, we can therefore choose 
qo = q. p will then be admissible if ip(r a , s) is a continuous function. 

J a was originally introduced as something with compact support (in r a ), so A, B a , and C must also have compact 
support (in 7'"). Using this along with the requirement that p vanish outside W implies that ip = g/47r|r| in this 
region. So any choice of ip which satisfies this 'boundary condition' and is continuous will generate a physically 
reasonable charge density. 

For the 3-current j a :— e™ J a , the (necessarily regular) portions contributed by <p can be temporarily ignored to 
show that 



N- 1 



B a + (...)- vPd, 



qr 



(A37) 



The divergence of this last term is equal to —qv a d a 6 3 (r) , which exactly cancels the divergence of B a given in l|A27() . 
It is therefore useful to write B a in the form 



B a =H a + v^da 



qr 
47r|r| ; 



(A38) 
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with H a an arbitrary piecewise continuous vector field satisfying d a H a — 0. This ensures that j a is nonsingular. 
Writing out the current in full, 

j a = N- 1 [H a + v^dp (r» - r a <p} , (A39) 
it is clear that j a will have support W if 

- (i^) < A40 » 

outside W. This also guarantees that supp(B Q ) = W, as required. 

This completes our study of the current moments. Essentially all physically interesting expansions can be extracted 
from a J a of the form l|A24|l . A, B a and C are all functions of (r a , s), and have support W. C has the form l|A30(l . 
where <p is an arbitrary continuous function equalling g/47r|r| 3 outside of W . (| A3 1(1 shows that A is also derived from 
ip (with qo — q). B a has the form (| A38|) . where H a is a piecewise continuous (3-) vector field satisfying d a H a = 0. 
The relevant portions of these results are also summarized in Sec. IIIII 



2. Stress-Energy Moments 

Let the multipole moments of T ab be denoted by the set {t bc , t abc , £ ai ■ Q ™ bc , . . .}, where t^ 1 - a ^ bc = t<*i-o»(M = 
i ai '" a " 6c . As with the current moments, it is convenient to keep track of this collection with a generating function 

00 <—i\ n 

G ab (k, s) := ^r fc d ' ' ' k c J Cl '" Cnab (s) . (A41) 

l|2.4ll inextricably links T ab and J a , so an expression for (T ab , (j> a b) as simple as the one for ( J a , cf> a } is not possible 
while retaining simple constraint and evolution equations. Still, one might expect that (T ab , (j> a b) should at least be 
proportional to fds (f ab , (j) ab ) (where f ab := (27r) 4 G ab ). Define a distribution $ afc = $( afc ) to make up the difference: 

{T a \ ^ ab ) = Jds (f ab (r, s) + $ afc (r, a), <f> ab (x)) . (A42) 

As before, we call all sets {t "} satisfying l)A41|l and (|A42p 'multipole moments of T^.' This is not a unique 
definition, however (even if $ were given ). Following [13,121, we will now pick out a set which very simply and 
naturally implies (|2.4|) . 

Using l|A14|) as a guide, (d a T ab , (frb) can be found by substituting a test function of the form </> Q t, = d a (pb into l|A42() . 
It is clear that the resulting expression depends on k a G ab , which is analogous to what happened when computing 
d a J a - In that case, (|A12|> showed that k a G a depended only on Q a . (|2.3|l is usually interpreted as an expression of 
global charge conservation, so it was not unreasonable that it only restricted the monopole moment. In the case of the 
stress-energy tensor, we expect that (|2.4H should have something to say about both the linear and angular momenta 
of the body (i.e. its monopole and dipole moments). We therefore suppose that k a G ab involves only t ab and t abc . This 
can be accomplished by letting 

f (oi-a„6)c = q ( A43 ) 



for all n > 2. 

Using this constraint implies that 



(d b T ab ^ a ) oc J ds (ik b G ab (k,s),<f> a (k)e~* k - z{s) ^ 

t ab d a fo(z( S ))+t abc d a d b c (z(s)) 



ds 



(A44) 
(A45) 



Dixon found that the first two moments can be given the special forms 20, 23] 

t ab = p (a v b) j ( A4 g) 
t abc = S a(b v c) ^ ( A47 ) 
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where we call p a the linear momentum, and S ab = S^ ab ^ the angular momentum. Using these expressions, 



(d b T ab , <j> a ) ex / ds 



(A48) 



In the absence of an electromagnetic field (or a current), the left-hand side of this equation vanishes, and the propor- 
tionality sign becomes an equality (since <& ab vanishes in this case) . Varying <j) a then recovers the standard equations 
of motion for a free particle: 



P a = 0, 
S ab = 2pt a v b] 



(A49) 
(A50) 



The situation is of course much more complicated when a field is present. Writing out ( x d a T ab , (fit) in full, l|2.4|) 
implies that 



p a 4> a + ~ (s ab - 2^1) d a 4> b + (ik b <S> ab ,4> a e- lk < s ^ \ = - (F afc J b ,0 a ) 



(A51) 



The right hand side of this equation can now be written in terms of the current moments. Using (|A10fl . the Fourier 
convolution theorem, and a Taylor series, it can be shown to be [2(J 



(F ab J b ,0 a ) = 
* Cl "- c " a (s) := 



1 



(2tt)* 

1 



E 

\p=0 



\n=0 

(-i) p 
p\ 



(-i) n 
n! 



/,/ ■■■l^Q dl - dpCl - 



b ,F a b {l)e- ll - z{s) \ 



(A52) 
(A53) 



Here, F ab does not quite represent the Fourier transform of F ab , which is not well-defined. It is instead equal to the 
Fourier transform of some function *F ab which coincides with the field in some neighborhood of S(s) n W, but has 
compact support. If *F ab is just as smooth as F ab , its precise form is irrelevant [2(i|. 

l|A52p now makes it natural to interpret the ^ "'s as multipole moments of the force density exerted on the body. 
We might therefore expect the net force to be proportional to and the net torque to $[ ob l. This identification can 
be made if <& ab has the form H3 



$ Qb (fc,s) = 



1 



4<( ab > - ik r 



(2tt)4 ^ 

Combining this with l|A51(l shows that 



^c(afe) _ l^J/ahc 

2 



E 

n=2 



(-0' 



km ' ' ' kr 



n -t- 2 

2vj/ c r-- c n( ab ) _ ~ iTi(ci---c n ab) 

n + 1 



f ds (p a + * Q ) cf> a + ~ (5 ab - 2p^ a v b] + 2^[ Qb l) d a (, 



= 



(A54) 



(A55) 



Given that this must hold for all possible choices of 4> a , it follows that 

p a = -* a , (A56) 
gab = 2 (p( a v b] ~ ¥ ab ^ . (A57) 

By construction, these are the only evolution equations implied by (|2.4|) . If we impose one more constraint equation: 

for n > 3, the chosen moments are unique in an appropriate sense. They are also sufficiently general to describe all 
physically interesting stress-energy tensors j2^|. Note that $ ab only depends on J a and F ab , and that the constraint 
equations (|A43p and l|A58|) are independent of these quantities. Portions of the stress-energy tensor which depend on 
the current have therefore been completely isolated from those which are not. 

Changes in the higher moments may once again be interpreted as 'equation of state' (this identification is actually 
more direct in this case). Their evolution is not completely arbitrary, however. Besides respecting the constraint 
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equations, they must also be chosen so that T ab remains physically reasonable. These extra restrictions would be 
provided by analogs of (|A36() and (|A39() . Although these will not be derived here, they can probably be constructed 
in a similar way. The presence of the field makes their derivation more complicated, but it should still be possible to 
repeat all of the steps carried out with the current moments. 

It suffices to note that for our purposes, this procedure has resulted in particularly natural definitions for the 
stress-energy moments - most importantly the linear and angular momenta. It can be shown that the choices made 
here imply that these momenta are given by (|3.1|) and i|3.2fl 20] . Interestingly, the net force and torque do not depend 
on G ab in any necessary way. They are apparently as independent of the details of the body's internal structure as 
possible. 

Before moving on, we can gain some insight into the force moments defined by i|A53[) . Only the first two of these 
are important here, and it is straightforward to show that 

* a = (j b (r,s),F ab (x)) , (A59) 
* [ah] = ij c (r,s),r [a F b ^ c (x)^ . (A60) 

Applying (|A24(1 . these expressions take the more explicit forms 



d 3 r <^ NJ b F ab + 



0_ 
~d~s 



d 3 r{NJ c r a e^F b]c + ( <p 



47r|r| 3 
<l 



e P b r p F at 



47r|r| 3 



e2r 7 v [a F b ^ c 



d_ 

ds 



^ ■ 



47r|r| 3 



e7r^r a e [ °F b ^ 



(A61) 
(A62) 



which were derived by a different method in Sec. IIIII 

Although i|A61(l and l|A62|) are exact, they are rather difficult to interpret. Their meaning is made considerably 
more transparent if the field can be expanded in a Taylor series inside S(s) n W. Then l|A59|l and (|A60|I together 
with IjAlJl show that 



1 

L 

V ab (-s) =s J2^Q dl " deaC ( d ^---9 de F b L )) 



(A63) 
(A64) 



When F ab is approximately constant throughout the charge, we recover the Lorentz force law, p a — — ^ a ~ —qF ab vt,. 
Unfortunately, these series are not useful when the field varies considerably over S(s) n W . And this is exactly what 
the self-field does. 
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